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ABSTRACT 


Finite  element  analysis  (FEA)  is  a  well-established  methodology  in  struetural  dynamics. 
However,  optimization  and/or  probabilistic  studies  can  be  prohibitively  expensive  because  they 
require  repeated  FEA  of  large  models.  Various  re-analysis  methods  have  been  proposed  in  order 
to  calculate  efficiently  the  dynamic  response  of  a  structure  after  a  baseline  design  has  been 
modified,  without  recalculating  the  new  response.  The  parametric  reduced-order  modeling 
(PROM)  and  the  combined  approximation  (CA)  methods  are  two  re-analysis  methods,  which  can 
handle  large  model  parameter  changes  in  a  relatively  efficient  manner.  Although  both  methods 
are  promising  by  themselves,  they  can  not  handle  large  EE  models  with  large  numbers  of  degrees 
of  freedom  (DOE)  (e.g.  100,000)  and  design  parameters  (e.g.  50),  which  are  common  in  practice. 
In  this  paper,  a  new  re-analysis  method  is  proposed  where  the  original  CA  method  is  modified  to 
further  improve  its  efficiency,  especially  for  problems  where  a  large  number  of  modes  is 
required.  The  modified  CA  (MCA)  method  is  then  integrated  with  the  PROM  approach,  in  order 
to  formulate  a  new  efficient  re-analysis  method  that  is  suitable  for  large  PE  models  with  many 
design  parameters  that  vary  in  a  wide  range.  A  simple  frame  structure  is  used  to  explain  all  steps 
of  the  proposed  method.  Also,  a  vibro-acoustic  analysis  of  a  realistic  vehicle  PE  model  is 
presented  to  demonstrate  the  efficiency  and  accuracy  of  the  new  method.  A  design  optimization 
study  is  also  performed  to  highlight  the  accuracy  and  efficiency  of  the  proposed  re-analysis 
method. 
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1.  INTRODUCTION 


Finite  element  analysis  (FEA)  is  a  well-established  numerical  simulation  method  for 
structural  dynamics.  It  serves  as  the  main  computational  tool  for  NVH  analysis  in  the  low- 
frequency  range.  As  a  result  of  developments  in  numerical  methods,  as  well  as  advances  in 
computer  software  and  hardware,  FEA  can  now  handle  much  more  complex  models  far  more 
efficiently  than  even  a  few  years  ago.  However,  the  demand  for  computational  capabilities 
increases  in  step  with,  or  even  beyond  the  pace  of  these  improvements.  Eor  example,  automotive 
companies  are  constructing  more  detailed  models  with  millions  of  degrees  of  freedom  (DOE)  to 
study  vibro-acoustic  problems  in  higher  frequency  ranges.  Although  these  tasks  can  be  performed 
with  EEA,  the  computational  costs  are  enormous,  even  for  high-end  workstations  with  the  most 
advanced  software. 

There  are  also  other  factors  that  increase  the  computational  costs.  When  design  changes  are 
involved,  the  EEA  analysis  must  be  repeated  many  times  in  order  to  obtain  the  optimum  design. 
Eurthermore  in  probabilistic  analysis  where  parameter  uncertainties  are  present,  the  EEA  analysis 
must  be  repeated  for  a  large  number  of  sample  points.  In  such  cases,  the  computational  cost  is 
even  higher,  if  not  prohibitive. 

Re-analysis  methods  are  intended  to  analyze  efficiently  structures  that  are  modified  due  to 
various  changes.  The  objective  is  to  evaluate  the  structural  response  for  such  changes  without 
solving  the  complete  set  of  modified  analysis  equations.  Several  reviews  have  been  published  on 
re-analysis  methods  [1-3].  These  methods  are  usually  based  on  local  and  global  approximations. 
Eocal  approximations  are  very  efficient  but  they  are  effective  only  for  small  structural  changes. 
Global  approximations  are  preferable  for  large  changes,  but  they  are  usually  computationally 
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expensive  especially  for  cases  with  many  design  parameters.  The  well-known  Rayleigh-Ritz  re¬ 
analysis  procedure  [4,  5]  belongs  to  the  category  of  local  approximation  methods.  The  mode 
shapes  of  a  nominal  design  can  be  used  to  form  a  Ritz  basis  for  predicting  the  response  in  a  small 
parametric  zone  around  the  nominal  design  point.  However,  it  is  incapable  of  capturing  relatively 
large  design  changes. 

A  parametric  reduced-order  modeling  (PROM)  method,  developed  by  Balmes  [6,  7], 
expands  on  the  Rayleigh-Ritz  method,  by  using  the  mode  shapes  from  a  few  selected  design 
points  to  predict  the  response  throughout  the  design  space.  The  PROM  method  still  belongs  to 
the  category  of  local  approximation  methods,  but  it  can  handle  relatively  larger  parameter 
changes  because  it  uses  multiple  design  points.  An  improved  component-based  PROM  method 
has  been  introduced  by  Zhang  et.  al.  [8,  9]  for  design  changes  at  the  component  level.  The 
method  was  implemented  in  NASTRAN  DMAP  (Direct  Matrix  Abstraction  Program),  and  its 
efficiency  and  accuracy  were  demonstrated  using  a  medium  size,  finite-element  vehicle  model. 

The  PROM  method  using  a  ‘parametric’  approach  (Section  2.2),  has  been  successfully 
applied  to  design  optimization  and  probabilistic  analysis  of  vehicle  structures.  However,  the 
‘parametric’  approach  is  only  applicable  to  problems  where  the  mass  and  stiffness  matrices  can 
be  approximated  by  a  polynomial  function  of  the  design  parameters  and  their  powers.  A  new 
‘parametric’  approach  using  Kriging  interpolation  [10]  has  been  recently  proposed  [11].  It 
improves  efficiency  by  interpolating  the  reduced  system  matrices  without  needing  to  assume  a 
polynomial  relationship  of  the  system  matrices  with  respect  to  the  design  parameters  as  in  [6,  7]. 

The  Combined  Approximations  (CA)  method  [12-14]  combines  the  strengths  of  both  local 
and  global  approximations,  being  accurate  even  for  large  design  changes.  The  CA  method  is  a 


UNCLASSIFIED;  Dist  A.  Approved  for  public  release 


4 


combination  of  binomial  series  (local)  approximations  (also  called  Neumann  expansion 
approximations)  and  reduced  basis  (global)  approximations.  Originally,  the  CA  method  was 
developed  for  linear  static  re-analysis.  It  has  also  been  extended  to  eigen-problem  re-analysis  [15- 
19].  Accurate  results  and  significant  computational  savings  have  been  reported.  A  comparative 
study  [20]  has  shown  that  the  method  is  suitable  for  large  design  changes. 

All  reported  studies  on  the  CA  method  [12-19],  have  used  relatively  simple  frame  or  truss 
systems  for  static  or  dynamics  analysis  with  a  relatively  small  number  of  DOF  and/or  modes.  For 
these  problems,  the  computational  efficiency  was  improved  by  a  factor  of  5  to  10.  Such  an 
improvement  is  beneficial  for  a  single  design  change  evaluation  or  even  for  gradient-based 
design  optimization,  where  only  a  limited  number  of  re-analyses  (e.g.  less  than  50)  is  performed. 
However,  the  computational  efficiency  of  the  CA  method  may  not  be  adequate  in  simulation- 
based  (e.g.  Monte-Carlo)  probabilistic  dynamic  analysis  of  large  finite-element  models,  where  re¬ 
analysis  must  be  performed  hundreds  or  thousands  of  times  in  order  to  estimate  accurately  the 
reliability  of  a  design.  The  main  reason  that  the  CA  method  is  still  expensive  for  such  problems 
is  the  considerable  cost  for  calculating  the  basis  vectors  used  to  approximate  the  response. 

A  large  number  of  modes  must  be  calculated  and  used  in  a  dynamic  analysis  of  a  large  finite- 
element  model  with  a  high  modal  density,  even  if  a  reduced-order  modeling  approach  is  used 
(Section  2.1).  In  such  a  case,  the  implicit  assumption  of  the  CA  method  that  the  cost  of  solving  a 
linear  system  is  dominated  by  the  cost  of  matrix  decomposition  is  no  longer  valid  (see  section 
3.1).  As  a  result,  the  computational  savings  from  using  the  CA  method  may  not  be  substantial. 
For  this  reason,  a  modified  combined  approximation  (MCA)  method  is  proposed  in  this  paper 
and  is  integrated  with  the  PROM  method  for  improving  accuracy  and  computational  efficiency. 
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The  MCA  method  approximates  efficiently  the  modes  of  a  number  of  designs.  These  modes  are 
subsequently  used  by  the  PROM  method  in  an  optimization  process.  The  computational  savings 
can  be  substantial  for  problems  with  a  large  number  of  design  parameters.  Examples  demonstrate 
the  benefits  of  the  proposed  re-analysis  method. 

The  paper  is  organized  as  follows.  The  existing  PROM  and  CA  methods  are  reviewed  in 
Section  2  after  a  brief  introduction  to  reduced  order  modeling  for  dynamic  systems.  Section  3 
describes  the  proposed  MCA  method  and  its  integration  with  gradient-based  optimization  and  the 
PROM  method.  Section  4  demonstrates  the  advantages  of  the  proposed  MCA  method  and  its 
integration  with  optimization  and  the  PROM  method,  using  a  simple  frame  structure  example 
and  a  complex  vehicle  model.  Finally,  conclusions  are  drawn  in  Section  5. 

2.  REVIEW  OF  EXISTING  METHODS 

2.1  Reduced  Order  Modeling  for  Dynamic  Response 

For  an  undamped  structure  with  stiffness  and  mass  matrices  K  ,  and  M  respectively,  under 
the  excitation  force  vector  F ,  the  equations  of  motion  (EOM)  for  frequency  response  are 

[K-ft;^M]d  =  F  (1) 

where  the  displacement  d  is  calculated  at  the  forcing  frequency  (O .  If  the  response  is  required  at 
multiple  frequencies,  the  repeated  direct  solution  of  Eq.  (1)  is  computationally  very  expensive 
and  therefore,  impractical  for  large  scale  finite-element  models. 

The  reduced  order  model  (ROM)  of  this  Section  is  a  subspace  projection  method.  Instead  of 
solving  the  original  response  equations,  it  is  assumed  that  the  solution  can  be  approximated  in  a 
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subspace  spanned  by  the  dominant  mode  shapes.  A  modal  response  approach  can  be  used  to 
calculate  the  response  more  efficiently.  A  set  of  eigen-frequencies  <»,  and  corresponding 
eigenvectors  (mode  shapes)  (p,  are  first  obtained.  Then,  the  displacement  d  is  approximated  in 
the  reduced  space  formed  by  the  first  n  dominant  modes  as 

d  =  W ,  (2) 

where  O  =  [tpj  <P2  •••  <PJ  is  the  modal  basis  and  U  is  the  vector  of  principal  coordinates  or 

modal  degrees  of  freedom  (DOF).  Thus,  the  EOM  of  Eq.  (1)  can  be  transformed  from  the 
original  physical  to  the  modal  degrees  of  freedom  as 

[6^K6-ft»^6^M6ju  =  6^F.  (3) 

Solving  for  the  modal  response  U  and  projecting  it  back  to  the  physical  coordinates,  the 
response  d  can  be  recovered.  If  is  the  maximum  excitation  frequency,  it  is  common  practice 
to  retain  the  mode  shapes  in  the  frequency  range  of  The  system  modes  in  the  high 

frequency  range  can  be  safely  truncated  with  minimal  loss  of  accuracy. 

Due  to  the  modal  truncation,  the  size  of  the  ROM  is  drastically  reduced,  compared  to  the 
original  model.  However,  the  size  increases  with  the  maximum  excitation  frequency.  An  added 
benefit  of  the  ROM  is  that  Eqs.  (3)  are  decoupled  because  of  the  orthogonality  of  the  mode 
shapes.  Therefore  these  equations  can  be  solved  separately,  reducing  further  the  overall 
computational  effort. 

2.2  Parametric  Reduced-Order  Modeling  Method 
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The  parametric  reduced-order  modeling  (PROM)  method  is  a  subspace  approximation 
method  developed  to  improve  the  efficiency  of  re-analysis.  It  assumes  that  the  dynamic  response 
of  a  new  design,  which  is  obtained  by  relatively  large  parameter  changes  from  the  baseline 
design,  can  be  approximated  in  the  subspace  spanned  by  the  mode  shapes  of  some  representative 
designs. 

The  quality  of  the  approximation  depends  on  the  selection  of  the  representative  designs. 
Balmes  et  al.  [6,  7]  suggested  a  procedure  for  obtaining  such  designs.  Consider  a  design  defined 
in  terms  of  m  design  parameters  p  such  as  plate  thicknesses,  spot  weld  diameters,  spot  weld 
pitch  of  a  car  body,  etc.  The  ranges  in  which  these  parameters  and  their  baseline  values  can  vary 
are  estimated  by  the  designer.  We  consider  (m  -I- 1)  designs  where  the  first  design  corresponds  to 
the  baseline  values,  Po,  of  the  design  parameters  and  the  remaining  designs  are  obtained  by 
perturbing  the  design  parameters  from  their  baseline  values  to  their  upper  limits,  one  at  a  time. 
Let  P;,/  =  l,---,m  be  the  vector  of  parameter  values  of  the  design.  The  (m-i-l)  designs  are 
represented  by  (m-l-l)  points  in  the  space  of  the  design  parameters;  the  first  design  corresponds 
to  the  origin,  and  the  remaining  m  designs  correspond  to  points  located  on  each  of  the  m  axes. 

In  the  PROM  approach  [7],  the  mode  shapes  of  a  new  design,  with  stiffness  and  mass 
matrices  and  ,  are  approximated  in  the  subspace  P  spanned  by  the  modes  of  the  above 

(m-l-l)  designs, 

p=l<i>,  «>,  ...  ©J,  (4) 

where  Og  is  the  modal  matrix  composed  of  the  dominant  mode  shapes  for  the  design  point  with 

^  .7 

all  parameters  at  their  lower  limits  (baseline  design),  and  O,  is  the  modal  matrix  for  the  i 
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design  point,  which  is  the  same  with  the  baseline  design  except  that  its  i'"  parameter  is  at  its 
upper  limit.  For  m  design  parameters,  (m  +  l)  eigenvalue  problems  must  be  solved  in  order  to 

form  matrix  P .  Therefore,  both  the  cost  of  obtaining  modal  matrices  O,  and  the  size  of  matrix 

P  increase  linearly  with  the  number  of  parameters.  Because  basis  P  captures  the  dynamic 
characteristics  in  each  dimension  of  the  parameter  space  the  response  at  a  new  design  point  in  the 

parameter  space  can  be  approximated  accurately  if  Pis  used  in  Eq.  (3)  instead  of  ^ . 

PROM  uses  the  following  algorithm  to  compute  the  mode  shapes  of  a  new  design: 

1)  Select  lower  and  upper  bounds  of  the  design  parameters. 

2)  Find  the  modes  of  the  baseline  design  and  the  designs  corresponding  to  the  m  points  in 
the  design  space,  and  form  subspace  basis  P . 

3)  Condense  matrices  and  as 

KR=P^KpP  MR=P^MpP,  (5) 

4)  Perform  an  eigen-analysis  on  the  reduced  matrices  to  calculate  eigenvector  0  . 

3)  Reconstruct  the  approximate  eigenvectors  using 

®^=P0.  (6) 

The  computational  cost  of  PROM  is  very  high  when  the  number  of  parameters  m  is  large 
because:  a)  m  eigen-analyses  must  be  performed  to  determine  subspace  basis  P  in  Eq.  (4),  b)  P 
is  usually  a  very  high  density  matrix,  and  c)  the  cost  of  forming  the  condensed  matrices  (by 
computing  the  triple  matrix  products  in  Eq.  (5))  is  high. 

Balmes  et  al  [7]  have  proposed  a  ‘parametric’  approach  to  reduce  the  cost,  for  the  special 
case  where  matrices  and  can  be  approximated  by  the  n''  order  polynomial  functions  of 
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m  material  parameters  p,  with  matrices  K,^  and  as  coefficients.  For  most  parameters  of 

material  properties  used  in  finite  element  modeling,  a  third  order  polynomial  is  sufficient.  The 
condensed  stiffness  and  mass  matrices  can  be  calculated  as  follows 

m  n  .  ,  m  n  . 

K«=ZE(P"K,P)p/  M.=XE(p"M„p)p/.  (7) 

i=l  j=0  i=l  j=0 

Because  K..  and  M.. are  constant  matrices,  the  condensations  P^K.P  and  P^M.P  are 
performed  only  once  and  then  they  are  used  to  form  matrices  and  in  Eq.  (7).  The  cost  of 
matrix  condensation  of  Eq.  (5)  is  largely  reduced  because  this  equation  involves  only  a  sum  of 
matrices  of  reduced  size  for  any  design.  However,  the  cost  of  the  parametric  approach  is  still  high 
when  the  number  of  parameters  m  is  high  because  (m  + 1)  eigen-analyses  must  be  performed  to 
form  basis  P. 

With  the  ‘parametric’  approach,  PROM  requires  first  to  form  the  subspace  basis  P  in  Eq. 
(4),  and  second  to  identify  matrices  and  ,  and  establish  the  relationship  of  Eq.  (7).  The 

cost  of  these  steps  is  considered  as  overhead.  The  remaining  cost  of  each  re-analysis  in  PROM  is 
minimal.  The  following  two  observations  can  be  made.  First,  the  overhead  cost  is  proportional  to 
the  number  of  design  parameters,  increasing  therefore,  the  efficiency  of  the  PROM  method  when 
the  number  of  re-analysis  is  significantly  larger  than  the  number  of  design  parameters.  This  is  the 
case  in  simulation-based  probabilistic  analysis.  Second,  the  efficiency  of  the  PROM  approach 
can  be  improved  by  reducing  the  above  overhead  cost. 

A  Kriging  method  is  introduced  in  [11]  to  reduce  the  cost  of  identifying  matrices  and 

My  in  Eq.  (7),  by  interpolating  the  reduced  system  matrices.  In  this  paper,  a  modified  combined 
approximations  (MCA)  method  is  proposed  that  can  drastically  reduce  the  computational  effort 
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in  forming  P  for  large,  finite-element  models  with  a  large  number  of  modes.  The  MCA  method 
is  described  in  Section  3  after  a  brief  introduction  of  the  CA  method  in  the  next  subsection.  We 
will  demonstrate  that  integration  of  the  MCA  and  PROM  methods  realizes  considerable 
computational  savings  with  a  small  loss  of  accuracy. 

2.3  The  Combined  Approximations  (CA)  Method 

The  PROM  method  requires  an  eigenvalue  analysis  for  multiple  designs  (comer  points)  to 
form  a  basis  for  approximating  the  eigenvectors  at  other  designs.  This  Section  provides  the 
fundamentals  of  the  combined  approximations  (CA)  method  [15-19]  which  does  not  require 
eigenvalue  analyses  for  multiple  designs.  Instead,  a  subspace  basis  is  formed  through  a  recursive 
process. 

If  Kq  and  Mq  are  the  stiffness  and  mass  matrices  of  the  original  (baseline)  design,  the  exact 
mode  shapes  obtained  by  solving  the  following  eigen-problem 

Ko^o  (8) 

We  want  to  approximate  the  mode  shapes  of  a  modified  design  (subscript  p)  with  stiffness  and 
mass  matrices 

K^=Ko+AK  M^=Mo+AM  (9) 

where  AK  and  AM  represent  large  perturbations.  The  CA  method  estimates  the  new 
eigenvalues  and  eigenvectors  without  performing  an  exact  eigenvalue  analysis. 

The  eigen-problem  for  the  modified  design  can  be  expressed  as 

Op=ApK~^MpOp-K~^AKOp,  (10) 
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Leading  to  the  following  recursive  equation 

^  P,j  =  Pj-l 

which  produces  a  sequence  of  approximations  of  the  mode  shapes  ^ pj,  j  =  1,2,..., s  .  The  CA 
method  uses  the  changes  ^  -  to  form  a  subspace  basis  to  approximate  the  modes 

of  the  new  design.  In  order  to  simplify  the  calculations,  pQ> p  j_i  in  Eq.  (10)  is 

replaced  with  XpK.'^'MpQ}^  and  Eq.  (11)  becomes 

(12) 

Eq.  (12)  shows  that  the  basis  vectors  satisfy  the  following  recursive  equation 

R^.  =-K-'AKR._,  j  =  2,...,s,  (13) 

where  the  first  basis  vector  is  assumed  to  be  Rj  =  Kq'M^Oq  . 

The  CA  method  forms  a  subspace  basis 

R  =  [R,  R,  ■■■  R,]  (14) 

where  s  is  usually  between  3  and  6  [13-17].  The  mode  shapes  of  the  new  (K^,  M^)  design  are 
then  approximated  in  the  subspace  spanned  by  R^  using  the  following  algorithm: 

1)  Condense  the  stiffness  and  mass  matrices 

K^=R"K^R  M^=R"M^R,  (15) 

2)  Solve  the  reduced  eigen-problem  to  calculate  eigenvector  matrix  0  . 

3)  Reconstruct  the  approximate  eigenvectors  of  the  new  design  as 

d^=R0.  (16) 
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The  eigenvalues  of  the  new  design  are  approximated  by  the  eigenvalues  of  the  reduced  eigen- 
problem. 

The  CA  method  has  three  main  advantages;  a)  it  only  requires  a  single  matrix  decomposition 
of  stiffness  matrix  Kq  in  Eq.  (13)  to  calculate  the  subspace  basis  R  ,  b)  it  is  accurate  because  the 

basis  is  updated  for  every  new  design,  and  c)  the  eigenvectors  of  a  new  design  are  efficiently 
approximated  by  Eq.  (16)  where  the  eigenvectors  0  correspond  to  a  much  smaller  reduced 
eigen-problem.  However  if  a  large  number  of  re-analyses  are  needed,  the  computational  cost  can 
be  high  because  a  new  basis  must  be  calculated  in  CA  for  every  re-analysis.  Note  that  steps  (15) 
and  (16)  are  similar  to  steps  (5)  and  (6)  of  PROM.  The  former  use  basis  R  and  the  latter  use 
basis  P. 

The  CA  method  is  more  efficient  than  PROM  for  design  problems  where  few  re-analyses  are 
required  for  two  reasons.  Eirst,  it  does  not  require  calculation  of  the  eigenvectors  ,  /  =  l,---,m , 
of  the  m  comer  design  points,  and  second  the  cost  of  matrix  condensation  of  Eq.  (15)  is  much 
lower  than  that  of  Eq.  (5),  because  the  size  (number  of  columns)  of  basis  R  is  not  proportional 
to  the  number  of  parameters  m  as  in  basis  P. 

As  discussed  in  Section  2.2,  the  PROM  approach  is  efficient  only  when  the  ‘parametric’ 
relationship  is  already  established  because  a  large  overhead  cost,  proportional  to  the  number  of 
design  parameters,  is  required.  On  the  contrary,  the  CA  method  does  not  require  such  an 
overhead  cost  because  the  re-analysis  cost  is  not  proportional  to  the  number  of  design 
parameters.  The  CA  method  is  therefore,  more  suitable  than  PROM,  if  the  number  of  re-analyses 
is  less  than  or  comparable  to  the  number  of  design  parameters. 
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The  CA  method  can  however,  become  expensive  when  many  re-analyses  are  needed  because 
a  new  basis  R  and  the  condensed  mass  and  stiffness  matrices  in  Eq.  (15)  must  be  calculated  for 
each  re-analysis.  Examples  where  many  analyses  are  needed  are  optimization  problems  in  which 
a  Genetic  Algorithm  is  employed  to  search  for  the  optimum,  and  probabilistic  analysis  problems 
in  which  Monte-Carlo  simulation  is  used.  The  PROM  method  can  be  more  suitable  for  these 
problems  because  the  subspace  basis  R  does  not  change  for  every  new  design  point. 

In  the  literature,  the  accuracy  and  efficiency  of  the  CA  method  has  been  mostly  tested  on 
problems  involving  structures  with  up  to  few  thousands  of  DOEs,  such  as  frames  or  trusses  [12- 
19].  Moreover  to  the  best  of  the  authors’  knowledge,  it  has  not  been  integrated  with  commercial 
finite  element  codes.  We  have  tested  the  CA  method  using  the  structural  dynamics  response  of  a 
medium  size  (65,000  DOE),  finite-element  model.  Due  to  its  high  modal  density,  there  were 
more  than  two  hundred  dominant  modes  in  the  frequency  range  of  zero  to  50  Hz.  It  was  observed 
that  the  computational  savings  of  the  CA  method,  using  the  recursive  process  of  Eq.  (13),  were 
not  substantial.  Eor  this  reason,  we  introduce  in  this  paper  a  modified  combined  approximations 
method  (MCA)  by  modifying  the  recursive  process  of  Eq.  (13).  The  new  method,  which  is 
described  in  Section  3.1,  is  much  more  efficient  than  the  original  CA  method  for  large  size 
models. 


3.  A  NEW  METHOD  FOR  EFFICIENT  RE-ANALYSIS 

3.1  Modified  Combined  Approximation  (MCA)  Method 

The  cost  of  calculating  the  subspace  basis  in  Eq.  (13)  consists  of  one  matrix  decomposition 
(DCMP)  and  one  forward-backward  substitution  (PBS).  The  DCMP  cost  is  only  related  to  the 
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size  and  density  of  the  symmetric  stiffness  matrix,  while  the  FBS  cost  depends  on  both  the  size 
and  density  of  the  stiffness  matrix  and  the  number  of  columns  of  ®  o . 

As  we  mentioned  in  Section  2.1,  the  higher  the  frequency  range  of  interest  the  more  modes 
are  needed  to  predict  the  structural  response  accurately.  In  such  a  case,  although  a  single  DCMP 
is  needed  in  Eq.  (13),  the  number  of  columns  in  may  increase  drastically,  thereby  increasing 

the  cost  of  the  repeated  FBS.  When  the  number  of  dominant  modes  becomes  very  large,  the  cost 
of  performing  the  calculations  in  Eq.  (13)  can  be  dominated  by  the  FBS  cost.  For  example,  the 
vehicle  model  of  Section  4.2  (Figure  5)  has  about  65,000  degrees  of  freedom,  and  1050  modes  in 
the  frequency  range  of  0-300  Hz.  If  has  only  one  mode,  the  cost  of  one  DCMP  is  1.1  seconds 
and  the  cost  of  one  FBS  is  less  than  0.1  seconds.  In  this  case,  the  total  cost  is  dominated  by  the 
DCMP,  and  the  CA  method  reduces  the  cost  of  one  re-analysis  considerably.  However,  if  has 

1050  modes  the  cost  of  FBS  increases  to  29  seconds,  dominating  the  cost  of  the  DCMP. 

The  CA  method  can  therefore,  improve  the  efficiency  only  when  the  number  of  retained 
modes  is  small.  Otherwise,  the  computational  savings  do  not  compensate  for  the  loss  of  accuracy 
from  using  Kq  (stiffness  matrix  of  baseline  design)  instead  of  (stiffness  matrix  of  new 

design).  For  this  reason,  a  modified  combined  approximations  (MCA)  method  is  proposed  in  this 
paper. 

The  MCA  method  uses  a  subspace  basis  T  whose  columns  are  constructed  using  the 
recursive  process 

T,  =K-;(MJ,_,)  i  = 

instead  of  that  of  Eq.  (13).  The  selection  of  the  appropriate  number  of  basis  vectors  5  is 
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discussed  later.  The  only  difference  between  Eqs.  (13)  and  (17)  is  that  matrix  is  inverted  in 


the  former  while  matrix  is  inverted  in  the  latter.  The  DCMP  of  must  be  repeated  for 

every  new  design.  However,  the  cost  of  the  repeated  DCMP  does  not  significantly  increase  the 
overall  cost  in  Eq.  (17),  because  the  latter  is  dominated  by  the  EBS  cost.  Note  that  the  iterative 
process  of  Eq.  (17)  provides  a  continuous  mode  shape  updating  of  the  new  design.  If  the  process 
converges  in  s  iterations,  the  mode  shapes  will  be  the  exact  mode  shapes  .  However,  Eq. 

(13)  does  not  have  the  same  property.  Eor  this  reason,  the  vectors  T,  provide  a  more  accurate 
approximation  of  the  exact  mode  shapes  than  the  R,  vectors  of  the  original  CA  method. 
This  is  an  important  advantage  of  MCA. 

Because  the  mode  shapes  T,.  in  Eq.  (17)  can  quickly  converge  to  the  exact  mode  shapes  , 
only  one  iteration  (i.e.  5=1)  may  be  needed,  resulting  in 

T  =  Ti.  (18) 

If  the  convergence  is  slow,  multiple  sets  of  updated  mode  shapes  can  be  used  so  that 

T  =  [Oo  T,  T,  -  Tj.  (19) 

Eor  better  accuracy,  the  above  basis  can  also  include  the  mode  shapes  of  tho  baseline  design. 
The  approximate  modes  T,  are  more  accurate  than  the  CA  vectors  R,  in  approximating  the 
exact  mode  shapes  ,  because  modes  T,  are  updated  for  each  new  design.  MCA  can  therefore, 

achieve  similar  accuracy  to  the  CA  method  using  fewer  modes.  As  we  will  demonstrate  later 
using  a  vehicle  model,  the  MCA  method  achieves  good  accuracy  with  only  1  or  2  basis  vectors, 
whereas  the  CA  method  requires  3  to  6  basis  vectors  [13-17]. 

In  summary,  the  proposed  MCA  method  involves  four  steps  in  calculating  the  approximate 
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eigenvectors  O  ^ : 

1)  Find  basis  T  using  Eq.  (18)  or  Eq.  (19). 

2)  Calculate  the  condensed  stiffness  and  mass  matrices  and 

K^=T^K^T  =T^M^T.  (20) 

3)  Solve  the  following  reduced  eigen-problem  to  calculate  the  eigenvalues  and  the 
projections  of  the  modes  in  the  reduced  space  spanned  by  T 

(K« -/IM«)0  =  O.  (21) 

4)  Reconstruct  the  approximate  eigenvectors  ®  ^  as 

d^=T0.  (22) 

it  should  be  noted  that  the  slightly  increased  cost  of  using  Eq.  (19)  instead  of  Eq.  (18)  is 
usually  smaller  than  the  realized  savings  in  steps  (20)  through  (22)  due  to  the  smaller  size  of  the 
reduced  basis  T .  The  bases  of  Eqs.  (18)  and  (19)  are  smaller  in  size  than  the  CA  basis  of  Eq. 
(14)  for  comparable  accuracy.  The  MCA  method  requires  therefore,  less  computational  effort  for 
steps  (20)  through  (22).  The  computational  savings  compensate  for  the  increased  cost  of  DCMP 
for  dynamic  re-analysis  of  large  finite-element  models  with  a  large  number  of  dominant  modes. 

All  mode  shapes  in  Eq.  (21)  must  be  calculated  simultaneously  in  order  to  ensure  that  the 
approximate  mode  shapes  are  orthogonal  with  respect  to  the  mass  and  stiffness  matrices. 

However,  the  cost  of  estimating  the  mode  shapes  ^^in  Eqs.  (20)  to  (22)  may  increase  quickly 

(quadratically)  with  the  number  of  modes,  and  as  a  result,  the  MCA  method  may  become  more 
expensive  than  a  direct  eigen-solution  when  the  number  of  dominant  modes  exceeds  a  certain 
limit. 


UNCLASSIFIED;  Dist  A.  Approved  for  public  release 


17 


One  way  to  circumvent  this  problem  is  to  divide  the  frequency  response  into  smaller 
frequency  bands  and  calculate  the  frequency  response  in  each  band  separately  instead  of  solving 
for  the  frequency  response  in  one  step.  The  modal  basis  T  in  Eq.  (19)  is  divided  into  k  groups  as 

T  =  [t'  •••  7^=]  (23) 

where, 

T‘  =  [a»;  t;  t'  •••  t;J.  (24) 

Each  group  T' contains  roughly  nik  original  modes  from  and  their  corresponding 

improved  modes.  The  eigenvectors  of  the  new  design  are  calculated  using  T'  instead  of  T  in  Eqs. 
(20)  to  (22).  The  process  is  repeated  k  times  using  a  modal  basis  that  is  1/fc  of  the  size  of  the 
original  modal  basis.  All  k  groups  of  eigenvectors  are  then  collected  together  to  calculate  the 
frequency  response  of  the  new  design.  As  demonstrated  in  the  numerical  examples  section,  this 
reduces  the  cost  considerably  with  minimal  loss  of  accuracy. 

3.2  Integration  of  the  MCA  Method  with  Optimization 

We  have  mentioned  that  the  MCA  method  provides  a  good  balance  between  accuracy  and 
efficiency  for  problems  that  require  a  moderate  number  of  re-analyses,  as  in  gradient-based 
optimization  problems.  Eor  problems  where  a  large  number  of  re-analyses  is  necessary,  such  as 
in  probabilistic  analysis  and  gradient-free  (e.g.  genetic  algorithms)  optimization,  a  combination 
of  the  MCA  and  PROM  methods  is  more  suitable  than  the  MCA  method. 

In  this  paper,  we  have  integrated  the  MCA  method  with  the  gradient-based  optimizer  in 
NASTRAN  v2001.  Eigure  1  shows  a  flowchart  of  the  optimization  process.  The  DMAP  (Direct 
Matrix  Abstraction  Program)  capabilities  in  NASTRAN  have  been  used  to  integrate  the  MCA 
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method  and  the  NASTRAN  modal  dynamic  response  and  optimization  (Sol.  200).  The 
highlighted  boxes  indicate  modifications  to  the  NASTRAN  optimizer.  Section  4.2  demonstrates 
how  this  process  was  used  to  optimize  the  vibro-acoustic  behavior  of  a  65,000  DOF,  finite- 
element  model  of  a  truck.  Using  the  MCA  method,  the  computational  cost  of  the  entire 
optimization  process  was  reduced  in  half  compared  with  the  existing  NASTRAN  approach. 


■^^riginal  Desig^^ 
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Sensitivity  Calculation 
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Improved  Design 
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Mode  Shapes  from  MCA 

I 
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I 

Recover  of  Response 

No 
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^hecjj 
Yes 
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Figure  1 .  Flowchart  for  MCA-enhanced  NASTRAN  optimization 


3.3  Integration  of  MCA  and  PROM  Methods 

The  PROM  method  requires  exact  calculation  of  the  mode  shapes  for  all  designs 
corresponding  to  the  corner  points  of  the  parameter  space  in  order  to  calculate  the  subspace  basis 
O .  The  required  computational  effort  can  be  prohibitive  if  a  large  number  of  parameters 
(optimization  design  variables)  is  used.  This  effort  can  be  reduced  substantially  if  the  modes  of 
each  corner  point  are  approximated  by  the  MCA  method.  In  this  case,  an  exact  eigen-solution  is 
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required  only  for  the  baseline  design. 


The  following  steps  describe  an  algorithm  to  integrate  the  MCA  and  PROM  methods: 

1.  Perform  an  exact  eigen-analysis  at  the  baseline  design  point  Pq  to  obtain  the  baseline 
mode  shapes, 

2.  Use  the  MCA  method  at  design  point  p, ,  where  all  parameters  are  set  at  their  low  limits 

except  the  parameter  which  is  set  at  its  upper  limit.  For  point  p^ ,  all  parameters  are  at  their 

lower  limits.  Obtain  approximate  mode  shapes  for  the  i'"  comer  point  using  the  following 
recursive  process 


=  7  =  2,3,-,/ 


(25) 


3.  Form  the  subspace  basis  T  as 

T  =  K  T„,  T,,  T,.,J  (26) 

where  m  is  the  total  number  of  parameters. 

4.  Obtain  the  approximate  mode  shapes  using  the  subspace  projection  procedure  of  Eqs. 

(5)  and  (6),  where  T  is  used  instead  of  P.  can  be  subsequently  used  in  a  modal  dynamic 

response  solution.  Only  step  4  is  repeated  in  re-analysis.  The  computational  savings  can  be 
substantial  especially  for  problems  where  many  re-analyses  are  needed,  such  as  in  probabilistic 
analysis  of  structural  dynamic  systems  using  Monte-Carlo  simulation. 


4.  NUMERICAL  EXAMPLES 


A  simple  frame  model  is  first  used  to  compare  the  accuracy  and  efficiency  of  the  MCA 


UNCLASSIFIED;  Dist  A.  Approved  for  public  release 


20 


method  with  that  of  the  original  CA  method.  Subsequently,  a  vehicle  model  is  used  to 
demonstrate  the  advantages  of  the  combined  MCA  and  PROM  method  in  optimizing  the  vibro- 
acoustic  response  of  a  vehicle.  Both  examples  are  performed  on  a  SUN  ULTRA  workstation 
using  NASTRAN  v2001.  The  CA,  MCA  and  PROM  methods  are  implemented  in  NASTRAN 
DMAP. 


4.1  A  Simple  Frame  Model 

The  MCA  method  is  applied  on  the  simple  eight-story  frame  structure  of  Figure  2.  The  same 
frame  model  has  been  extensively  used  in  the  literature  on  the  CA  method  [16-18].  A  lumped 
mass  representation  is  used  with  each  story  girder  having  a  mass  M.  The  girders  are  assumed  to 
be  non-deformable  with  the  initial  lateral  stiffness  {eI  /  )  of  50.0  (kips/inch)  for  stories  1  and 

2,  55.0  for  stories  3  to  6  and  60.0  for  stories  7  and  8.  The  corresponding  girder  lengths  are  2.3208 
ft  for  stories  1  and  2,  2.2482  ft  for  stories  3  to  6  and  2.1840  ft  for  stories  7  and  8. 

-  1 
-  2 
~  3 

-  J 
5 

►  6 

-  7 
^  8 

Figure  2.  Eight- story  frame  model 
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In  [16-18],  the  Young’s  moduli  of  all  girders  are  treated  as  design  parameters,  assuming 


implieitly  that  among  different  designs,  the  stiffness  matrix  varies  linearly  with  the  Young’s 
moduli  and  the  mass  matrix  remains  unehanged.  In  this  example,  we  allow  the  girder  lengths  to 
vary  resulting  in  a  nonlinear  ehange  of  the  stiffness  matrix  and  a  linear  ehange  of  the  mass  matrix 
from  the  baseline.  This  inereases  the  complexity  of  the  re-analysis. 

For  the  modified  design,  each  girder  mass  is  increased  by  25  percent,  and  each  girder 
stiffness  is  increased  by  a  factor  of  1.953.  Table  1  shows  the  first  mode  shapes  (p^  and  (p^  for 

the  modified  and  baseline  designs  respectively  obtained  from  NASTRAN.  These  results  are 
considered  exact  in  this  study. 


Table  1.  Accuracy  of  mode  shapes  for  the  frame  model 


Story 

# 

9,, 

t, 

^2 

Disp. 

Err(%) 

Disp. 

Err(%) 

Disp. 

Err(%) 

Disp. 

Err(%) 

1 

1.0000 

1.0000 

0 

1.0000 

0 

1.0000 

0 

1.0000 

0 

2 

0.9607 

0.9700 

0.97 

0.9592 

0.16 

0.9603 

0.04 

0.9583 

0.25 

3 

0.8836 

0.9109 

3.09 

0.8787 

0.55 

0.8825 

0.12 

0.8761 

0.85 

4 

0.8235 

0.8199 

0.43 

0.8153 

0.99 

0.8220 

0.18 

0.8149 

1.03 

5 

0.7407 

0.6954 

6.12 

0.7286 

1.64 

0.7388 

0.26 

0.7313 

1.27 

6 

0.6376 

0.5424 

14.93 

0.6219 

2.46 

0.6355 

0.34 

0.6285 

1.43 

7 

0.5170 

0.3673 

28.97 

0.4999 

3.32 

0.5149 

0.42 

0.5109 

1.20 

8 

0.2673 

0.1883 

29.53 

0.2563 

4.09 

0.2659 

0.51 

0.2620 

1.98 

The  MCA  method  was  used  to  approximate  the  mode  shapes  of  the  modified  design,  using 
the  recursive  process  of  Eq.  (23)  with  one  or  two  iteration  steps  (5'  =  lor  2).  In  Table  1,  tj  and 
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t2  are  the  once  (5  =  1)  and  twice  (5  =  2)  updated  mode  shapes.  Mode  shape  t2  is  very  close  to 


the  exact  mode  shape  (p^.  The  maximum  error  between  the  two  is  only  0.5%.  A  subspace  basis 
T'  =  Wo  t,]  (similar  to  Eq.  19)  was  also  tested.  The  resulting  approximate  mode  shape  is 
indicated  by  in  Table  1.  While  the  accuracy  of  tp^  is  not  as  good  as  1 2,  it  is  more  accurate 
than  tj,  indicating  that  the  basis  of  Eq.  (19)  provides  better  accuracy  than  that  of  Eq.  (18)  as 
expected. 

Eigure  3,  compares  tj  (MCA  with  once-updated  mode  shapes),  CA  with  twice-updated 
mode  shapes,  and  (p^  (exact  analysis  using  NASTRAN).  The  baseline  mode  is  also  shown.  It  is 
observed  that  tj  is  closer  to  (p^,  than  the  twice-updated  mode  from  the  CA  method.  This  supports 
the  claim  of  Section  3.1  that  the  MCA  basis  T,  provides  more  accurate  results  than  the  CA  basis 

R, 


Eigure  3.  Mode  shape  comparison  among  NASTRAN,  MCA  and  CA  methods 
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4.2  A  Vehicle  Model 


Figure  4  shows  a  65,000  DOF  finite-element  model  of  a  pickup  truck.  The  sound  pressure 
level  at  the  driver’s  ear  location  is  calculated  using  a  vibro-acoustic  analysis.  The  structural 
forced  vibration  response  due  to  unit  harmonic  forces  in  x,  y,  and  z  directions  at  the  engine 
mount  locations,  is  coupled  with  an  interior  acoustic  analysis.  The  first  and  second  eigen- 
frequecies  of  the  acoustic  volume  inside  the  cabin  are  95.9  Hz  and  128.3  Hz.  The  sound  pressure 
level  is  calculated  in  the  80  to  140  Hz  frequency  range. 

For  this  vibro-acoustic  analysis,  the  structure  and  fluid  domains  are  coupled  through 
boundary  conditions  ensuring  continuity  of  vibratory  displacement  and  acoustic  pressure.  A 
finite-element  formulation  of  the  coupled  undamped  problem  yields  the  following  system 
equations  of  motion  [21] 


where  the  primary  variables  are  the  vibratory  displacement  and  the  acoustic  pressure  p^.  The 
finite-element  representation  of  the  two  domains  consists  of  stiffness  and  mass  matrix  pairs 
(K  j  ,M5  )  and  ,M^ ) ,  respectively.  The  air  density  and  wave  speed  are  and  Cq  .  External 

forces  are  shown  in  the  right  hand  side  of  Eq.  (27). 

The  spatial  coupling  matrix  indicates  coupling  between  the  two  domains.  Due  to  this 
coupling,  which  is  usually  referred  to  as  “two-way  coupling”,  the  combined  structural-acoustic 
system  of  equations  is  not  symmetric.  If  the  acoustic  effect  on  the  structural  response  is  small,  the 
coupling  term  can  be  omitted,  resulting  in  the  so-called  “one-way  coupling,”  where  the  structural 
response  is  first  calculated  and  then  used  as  input  (f^  in  Eq.  27)  to  solve  for  the  acoustic 
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response. 


The  eoupled  strueture-aeoustie  system  ean  be  solved  either  by  a  direet  method,  or  more 
effieiently  by  a  modal  response  method  as  deseribed  in  Seetion  2.1.  A  modal  response  method 
ean  be  applied  to  both  the  struetural  and  the  aeoustie  domains. 


Figure  4.  A  piekup  truek  model 

4.2.1.  Combined  MCA  and  PROM  Methods 

A  re-analysis  was  performed  for  a  modified  design  for  whieh  the  following  five  plate 
thiekness  parameters  vary;  ehassis  and  its  eross  links,  eabin,  truek  bed,  left  door,  and  right  door. 
All  parameters  were  inereased  by  100%  of  their  baseline  values. 

The  sound  pressure  at  the  driver’s  ear  was  ealeulated  using  “two-way”  eoupling.  A  struetural 
modal  frequeney  response  analysis  was  used.  The  aeoustie  response  was  ealeulated  using  a  direet 
method  beeause  the  size  of  the  aeoustie  model  is  relatively  small.  For  the  struetural  analysis, 
1050  modes  were  retained  in  the  0  to  300  Hz  frequeney  range.  The  eombined  MCA  and  PROM 
approaeh  was  eompared  against  the  NASTRAN  direet  solution  for  a  modified  design  where  all 
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five  parameters  were  at  their  upper  limits.  Only  one  iteration  was  performed  in  Eq.  (17)  in  order 
to  get  one  set  of  once-updated  mode  shapes  for  each  corner  design  point.  The  subspace  basis, 
which  includes  information  for  all  five  design  parameters,  is  therefore,  represented  by 

T=K  Tji  T,,  ■■■  T,,].  (28) 

Figure  5  shows  the  percent  error  in  natural  frequencies  as  predicted  by  the  combined  MCA 
and  PROM  method  and  NASTRAN.  The  maximum  error  is  less  than  0.45%  in  the  entire 
frequency  range.  Figure  6  indicates  that  the  sound  pressures  calculated  by  both  methods  are 
almost  identical.  The  computational  effort  for  the  MCA  method  to  obtain  approximate  mode 
shapes  at  each  comer  design  point  is  about  30  seconds.  In  contrast,  it  takes  about  180  seconds  for 
an  exact  eigen- solution  using  NASTRAN. 

The  computational  costs  of  constructing  the  reduced  basis  (P  in  PROM  method  and  T  in 
PROM-i-CA  method)  are  compared  in  Table  2.  For  the  case  of  5  design  parameters,  the  total  cost 
was  reduced  from  1080  seconds  to  330  seconds.  The  computational  saving  is  more  significant  if 
the  number  of  design  parameter  increases. 
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Figure  5.  Percent  error  of  natural  frequencies  between  combined  MCA  and  PROM  method  and 


NASTRAN 


Figure  6.  Comparison  of  sound  pressure  at  driver’s  ear  between  combined  MCA  and  PROM 

method  and  NASTRAN 


Table  2.  CPU  time  for  constructing  reduced  basis 


Method 

Solving  for  mode 
shape  Oq  at 
baseline  design 

Solving  for  mode 
shapes  at  5 
corner  design 
points 

Total  Cost 

PROM 

180  sec. 

180*5=900  sec. 

1080  sec. 

PROM+MCA 

30  sec 

30*5=150  sec 

330  sec. 

4.2.2  Optimization  using  MCA  Method 

The  objective  here  is  to  minimize  the  sound  pressure  at  the  driver’s  ear.  In  order  to 
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demonstrate  the  capability  of  the  MCA  method  to  handle  a  large  number  of  parameters,  a  total  of 
41  design  parameters  are  used  representing  the  thickness  of  all  vehicle  components  modeled  with 
plate  elements.  All  thicknesses  are  allowed  to  change  by  100%  from  their  baseline  values.  Table 
3  gives  a  description  of  all  design  parameters.  All  design  parameters  are  at  their  low  bound  for 
the  optimization  initial  point. 


Table  3:  Description  of  design  parameters 


Prm. 

# 

Description 

(thickness  of) 

Prm. 

# 

Description 

(thickness  of) 

Prm. 

# 

Description 

(thickness  of) 

1 

Bumper 

15 

Radiator  mtg. 

29 

Tire,  front  right 

2 

Rails 

16 

Radiator  mtg.,  mid. 

30 

Tire,  rear  left 

3 

A-arm,  low  left 

17 

Ean  cover,  low 

31 

Tire,  rear  right 

4 

A- arm,  low  right 

18 

Ean  cover,  up 

32 

Engine  outer 

5 

A-arm,  up  left 

19 

Cabin 

33 

A-arm  conn.,  up  left 

6 

A-arm,  up  right 

20 

Cabin  mtg.  reinf. 

34 

A-arm  conn.,  up  right 

7 

Tire  rim 

21 

Door,  left 

35 

A-arm  conn.,  low  left 

8 

Engine  Oil-box 

22 

Door,  right 

36 

A-arm  conn.,  low  right 

9 

Ean 

23 

Bed 

37 

Glass,  left 

10 

Hood 

24 

Brake,  front  left 

38 

Glass,  right 

11 

Eender,  left 

25 

Brake,  front  right 

39 

Glass,  rear 

12 

Eender,  right 

26 

Rail  conn.,  rear 

40 

Glass,  front 

13 

Wheel  house,  left 

27 

Rail  mount 

41 

Rail  conn.,  front 

14 

Wheel  house,  right 

28 

Tire,  front  left 

The  MCA  method  approximates  the  mode  shapes  at  intermediate  design  points,  using  only 
Tj  in  Eq.  (17).  Therefore,  the  subspace  basis  at  each  optimization  step  is 

T  =  [0„  Tj.  (29) 
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Because  1050  modes  exist  in  the  frequency  range  of  0  to  300  Hz  of  the  initial  design,  the  size  of 


the  MCA  modal  basis  is  2100. 

The  cost  of  solving  for  1050  modes  directly  from  NASTRAN  is  180  seconds.  In  the  MCA 
method,  cost  of  solving  the  linear  system  of  equations  in  Eq.  (17)  is  31  seconds,  and  the 
additional  combined  cost  of  Eqs.  (20)  to  (22)  is  373  seconds,  resulting  in  a  total  cost  of  404 
seconds.  Eor  this  reason,  the  1050  modes  are  divided  into  21  groups  and  the  modes  in  each  group 
are  obtained  separately  as  explained  in  the  latter  part  of  Section  3.L  This  reduces  the  cost  of  Eqs. 
(20)  to  (22)  to  66  seconds.  The  total  cost  of  the  MCA  method  is  then  97  seconds,  which  is  about 
half  of  the  cost  of  the  direct  NASTRAN  method.  Details  of  the  computational  costs  are  provided 
in  Table  3. 


Table  3.  CPU  time  of  the  MCA  method 


k=l 

k=21 

Eq.  (17) 

31  sec. 

31  sec. 

Eq.  (20) 

258  sec. 

50  sec. 

Eq.  (21) 

48  sec. 

6  sec. 

Eq.  (22) 

67  sec. 

10  sec. 

Total  Cost 

404  sec. 

97  sec. 

The  optimization  process  of  Eigure  1  is  implemented  in  NASTRAN,  using  DMAP.  The 
gradient-based  optimizer  in  NASTRAN  (Sol  200)  needed  three  iterations  to  calculate  the  optimal 
design.  Eigure  7  compares  the  sound  pressure  at  the  driver’s  ear  between  the  optimal  design  and 
the  initial  design.  Eigure  8  indicates  the  percentage  increase  of  optimal  values  relative  to  the 
initial  values  for  all  41  design  parameters. 

The  Sequential  Quadratic  Programming  (SQP)  algorithm  of  NASTRAN  can  find  only  a 
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local  optimum  instead  of  the  global  optimum.  In  the  frequeney  range  of  80-140Hz,  the  maximum 
sound  pressure  is  slightly  redueed  from  7.9E-7  to  7.2E-7  Paseal.  Most  parameters  are  only 
slightly  ehanged.  The  largest  inerease  is  20%  for  the  rail  mount  thiekness  (parameter  #27). 


Eigure  7.  Comparison  of  sound  pressure  at  driver’s  ear,  between  initial  and  optimal  designs 


Eigure  8.  Pereent  inerease  of  optimal  design  parameters  relative  to  baseline  design  parameters 
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To  obtain  a  more  significant  design  improvement,  two  additional  studies  are  performed 
using  the  MCA  method  in  a  Genetic  Algorithm.  In  the  first  study,  20  initial  populations  and  4 
generations  were  used,  and  in  the  second  study  100  initial  populations  and  6  generations  were 
used.  Figures  9  and  10  show  that  the  optimization  results  are  affected  by  the  number  of  initial 
populations  and  the  number  of  generations.  While  more  number  of  initial  populations  and 
generation  produces  a  slightly  better  result,  both  studies  have  achieved  much  better  results  than 
the  SQP  algorithm.  In  the  case  of  100  initial  populations  and  6  generations,  the  sound  pressure  is 
reduced  from  7.9E-7  Pascal  to  2.0E-7  Pascal,  which  is  equivalent  to  about  15  dB  in  sound 
pressure  level  (SPL). 


Figure  9.  Comparison  of  sound  pressure  at  driver’s  ear,  between  baseline  and  optimal  designs, 

with  20  initial  populations  and  4  generations 
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Figure  10.  Comparison  of  sound  pressure  at  driver’s  ear,  between  baseline  and  optimal  designs, 

with  100  initial  populations  and  6  generations 


Figure  11.  Comparison  of  sound  pressure  at  driver’s  ear,  between  direet  NASTRAN  and  MCA 
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To  verify  the  accuracy  of  the  MCA  approximation,  the  sound  pressure  response  of  the 
optimal  design  from  MCA+GA  with  100  initial  populations  and  6  generations  was  evaluated  by 
both  direct  NASTRAN  and  MCA.  According  to  Figure  11,  the  MCA  method  is  very  accurate. 

For  a  comparable  to  MCA  accuracy,  the  original  CA  method  needed  three  sets  of  mode 
shapes  to  form  the  subspace  basis,  requiring  90  seconds  to  solve  the  linear  equations.  The  much 
larger  mode  basis  R  in  CA  increases  the  computational  cost  in  calculating  the  triple  matrix 
products  of  Eq.  (15).  Therefore  for  large  scale,  finite-element  models  with  a  high  modal  density, 
the  proposed  MCA  method  can  be  more  efficient  compared  to  a  complete  NASTRAN  analysis 
and  the  original  CA  method. 


5.  SUMMARY  AND  CONCLUSIONS 

An  accurate  and  efficient  modified  combined  approximations  (MCA)  re-analysis  method  has 
been  proposed  for  approximate  eigenvalue  analysis  of  large  size  models.  The  method  was  used  to 
improve  the  computational  efficiency  of  the  parametric  reduced-order  modeling  (PROM) 
method.  The  proposed  method  has  been  implemented  in  NASTRAN  and  demonstrated  on 
analysis  and  design  optimization  of  two  structures;  a  simple  frame  representing  a  building  and  a 
real-world  large-scale  model  of  a  truck.  We  showed  that  the  results  from  the  proposed  MCA 
method  are  very  similar  to  those  from  a  complete  NASTRAN  analysis  but  at  less  than  50%  of  the 
NASTRAN  cost. 
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